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Many-body effects on superfluidity and transition temperatures are calculated for optical lattices 
and uniform systems with ultracold multi-component Fermi gases. The induced interactions depend 
sensitively on the interactions between the multi-components and their densities. The s- and d- 
wave pairing gaps and critical temperatures are calculated for optical lattices at low (dilute) filling 
as well as near half filling to leading orders in the interaction strength with and without induced 
interactions. They can deviate strongly from dilute two-component systems and affect the phase 
diagram. In two dimensional optical lattices the induced interactions are singular at half filling and 
strongly affect s- and d-wave superfluidity. 
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Bardeen, Cooper, & SchriefTer's theory of supercon- 
ductivity (BCS) was early on extended with the many- 
body effect of induced interactions by Gorkov & Melik- 
Bharkudarov [H] who found that the BCS pairing gap 
and critical temperature was reduced by a factor ^e) 1 / 3 
in uniform three dimensional (3D) systems. Reductions 
have been observed in superconductors, 3 He Q, and in 
ultra-cold Fermi atom systems in traps and optical lat- 
tices 0, 0, H, 0| . The induced interactions were further 
discussed and extended to any number of spin states as 
well as including bosons [7]. Recently, the induced in- 
teractions were calculated numerically for lattices with 
two-component fermions within the Hubbard model, and 
near half filling they gaps and critical temperatures were 
strongly suppressed |8j]. Also induced interactions were 
calculated for uniform 3D 6 Li gases with three spin states 
1. 

It has become experimentally feasible to vary the inter- 
actions strength between states in atomic traps and op- 
tical lattices by using Feshbach resonances and to study 
the crossover from BCS to molecular BEC in mixtures of 
Fermi (and Bose) gases. Recently, multi-component sys- 
tems of Fermi atoms with several atomic states of inter- 
nal hyperspin have been produced such as 6 Li [10], 137 Yb 
with six nuclear spin states and heteronuclear mix- 
tures of 40 K and 6 Li [12]. Such multi-component systems 
have intriguing similarities with neutron, nuclear and 
quark matter where, for example, color superconductiv- 
ity between the 2 spin, 8 color and 2-3 flavor states may 
occur [13]. Ultracold Fermi atoms provide a highly con- 
trolled environment in which superfluidity can be studied 
in detail over a wide range of multi-component densities, 
temperatures, interaction strengths, etc. Additionally, 
2D optical lattices have the important quality that by 
changing the on-site interaction from attractive to repul- 
sive the BCS s-wave pairing changes to d-wave pairing, 
which is believed to be responsible for high temperature 
superconductivity 2D optical lattices may therefore 
also emulate high temperature superfluidity as well as 
Mott insulators, anti-ferromagnetic and other interesting 
phases and transitions (Tol. [la. [l~7T| . 



The purpose of this work is to calculate the s- and 
d-wave pairing gaps and critical temperatures for multi- 
component Fermi gases in optical lattices and uniform 
gases with and without induced interactions. For this 
purpose we employ the Hubbard Hamiltonian on a D- 
dimensional lattice, which describes optical lattices that 
are sufficiently deep for a one-band model to apply, 

H = U at<T m ia h il7 ' - t a\ a a, ja . (1) 

i,a<a' (ij),<? 

Here a\ a is the Fermi creation operator of the hyperspin 
states er = 1,2, ...,v, ni a = a\ a ai a the density and (ij) 
denotes nearest neighbours with hopping parameter t. 

We shall investigate s-wave pairing between to spin 
states, e.g. 1 and 2, due to an attractive on-site in- 
teraction Ui_2 < in the presence of multi-components 
j = 3,...,v with onsite interactions Uij. The mean field 
gap equation for singlet superfluidity at zero temperature 

is EE! 

A p - = 4E^(p',p)A p ^, (2) 
p p 

where M is the number of lattice points, E p = 
^J(e p - m) 2 + A p with e p = 2^ j= i,_d( 1 ~ c °sPi), and 
the chemical potentials are assumed equal /ii = /i2 
for the pairing components. The effective interaction 
Uil^ = Uu + Ui n d(p',p) includes induced interactions 
which will be calculated below. The density n = 1 — 
J2 p (e P /E p )ta,nh(E p /2k B T)/M, is also the filling frac- 
tion in units where the lattice constant is unity. Momenta 
are in the first Brillouin zone only \pi\ < it. 

Due to particle-hole symmetry results also apply re- 
placing the density by (2 — n) and chemical potentials 
by (4Dt — /i) , where ADt is the bandwidth in D dimen- 
sions. At low filling e p = tp 2 and results also apply to 
uniform systems replacing t = l/2m, where m is the par- 
ticle mass. 

We first calculate the zero temperature gaps for weak 
attraction |C/ 12 | <C t excluding induced interactions. The 
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s-wave gap A p / = Ao is then momentum-independent 
because the interaction Ui 2 is. The gap equation 
reduces to: 1 = — (Uq/M) Y^, p 1/2-Ep, from which the 
gap can be calculated analytically at low and half filling 
to leading orders in 3D and 2D (for ID see, e.g. (l9|). 

The 3D lattice has a critical coupli ng U c = 
-M/(£ q l/2e q ) = -8V2t/[J2Z P2i(VW8)(2l - 
1)!!/2 2 'Z!] ~ -7.913*, Q where a two-body bound state 
can be formed at zero density (fi\ — 0). This thresh- 
old corresponds to the unitary limit of infinite scattering 
length in the uniform system (see e.g. [21(), and natu- 
rally enters the sum in the gap equation. The dilute gap 
becomes to leading orders 



A 3 D = — m exp 
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n < 1, 
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at low densities where fi\ = tk F1 -C t. The level den- 
sity at the Fermi surface N(/j,i) = kF,i/(Air 2 t) enters in 
the exponent as in standard BCS theory. U c acts as a 
cutoff in the one-band Hubbard model which in the uni- 
form (continuum) limit is absorbed into the scattering 
lengths ajj between states i and j. Thus Eq. ((3J) is the 
finite lattice equivalent of the gap in the uniform system 
A 3 D = (8/e 2 ) A i 1 exp(7r/2fc Fil ai 2 ), with U 12 U C /(U C -U 12 ) 
replaced by Airai 2 /m such that threshold £7jj = U c cor- 
responds to \dij\ = 00. 

Near half filling the level density is surprisingly con- 
stant N(fii) ~ 0.143/i in a wide range At < fii < 8t 
around half filling. The gap can therefore is to leading 
orders 



Aq D = erf exp 



1 



N(jii)U : 
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where the prefactor a ~ 6.544.. is calculated numerically. 

The 2D lattice has a superfluid s-wave gap that be 
calculated from the gap equation [22] 



A 2 / 3 = y/8Tvmt exp 



Ant 

Th2 



n < 1, 



(5) 



when the density n = k F /2ir is small. The 2D level den- 
sity is N(fii) = l/(Ant). Eq. ((SJ assumes that the on-site 
coupling and the gap are small such that Ag^^/xi. This 
is not fulfilled at sufficiently low densities where instead 
Aq D = ATrtexp(8nt/Ui2), which also is the two-body 
binding energy in 2D at zero density. There is always 
a two-body bound state in 2D with purely attractive in- 
teraction and therefore U c vanishes. In 3D a similar pair 
condensate (a molecular BEC) appears when ai2 > 
corresponding to U\ 2 < U c . 

For intermediate fillings the pairing is more com- 
plicated. Near half filling n ~ 1 the level density 
N(e) = ln(16t/|e- At\)/(2n 2 t) has a logarithmic singu- 
larity due to the van Hove singularity. Calculating the 
r.h.s. of the gap equation therefore gives a double log: 
1 = |C/ 12 |/(47r 2 )ln 2 (32t/AF), to leading logarithmic or- 
ders, resulting in the gap [22] 



A 2 D = 32texp(-2^^^/t/\^h 2 ~\) , n ~ 1. 



(6) 



Induced interactions were originally calculated for 
dilute 3D two-component systems by Gorkov & Melik- 
Bharkudarov [l| and extended to multi-component Fermi 
gases also with Bose atoms • The induced interactions 
arise from two types of second order diagrams 0, 0, 0| 
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/ is the Fermi distribution of £j(q) = e q — [ij, k' = 
p + p' and k = p — p', where p.; are the momenta of 
the two pairing spin states. Induced interactions due to 
particle-hole loop diagrams from j = 3, v spin states 
are responsible for the second sum in Eq. ([7]) and has 
opposite sign. Near half filling the induced interactions 
are enhanced considerably because Umklapp processes 
enhance the phase space around where the denominator 
of Eq. |(7j) is singular. 

For later convenience we write the induced interactions 



U tnd = U^N(jjh)x(ui) - J2 U U U 2J N (^)x(l 

3=3 
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in terms of a screening function x of wi = k'/2fci? 1 and 
ujj = k/2kp,j- At low density the screening function 
reduces in 3D to the Lindhard function x 3Z5 (cj) = h + 



hi 



1+LU 

1-u 



whereas in 2D 



x 2D M = 



1- y/T^y: 



UJ < 1 
(J > 1 



(9) 



For weak interactions the gap is small and pairing oc- 
curs near the Fermi surface only. To calculate the s- 
wave gap we can therefore replace x(w) in Eq. ([7]) by 
its momentum average x(Mi> /^j) over the Fermi surface 
l e pl = l e p'l = Mi [H3]- At low densities pbj — tkpj and 
uoj = (kF,i/kF,j)s/(l + cos0)/2, where 8 is the angle be- 
tween p and p'. Averaging over this angle in Eq. J9]) 
gives 



X 2Z Vi,Mj) = 



1 



k pj j k 
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The screening function x 2D increases with multi- 
component filling only up to equal filling where a plateau 
appears. The result X 2I) (^i, H\) = 1 was earlier found by 
Nishida & Son [24] for uniform systems. In 3D uniform 

systems \ 3D = 1 when fj,j > fa, x 3D = (l/3)log(4e) 
when fij = fix and x 3D = (Mi/^Mi) log(/ii//Xj) when 

fij < a*i- Hi- 

X 2D (Mi) Hj) is shown in Fig. 1 for general filling. The 
plateau x 2D — 1 extends almost to half filling. Due to 
particle-hole symmetry the screening function is symmet- 
ric when any chemical potential [i is replaced by ADt — /1. 
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Figure 1: Averaged screening function X* D ((ii, flj) 
in a 2D optical lattice vs. multi-component chemical 
potentials. 



We find that the induced interactions has three loga- 
rithmic singularities at fii = 0,At,8t all at — At. 
The screening function itself is finite at these points. 
The singularity at half fillings was calculated numeri- 
cally by Kim et al. [8J] for a two component system. 
It can be calculated analytically by observing that near 
half filling the Fermi surface of components 1 and 2 is 
a square in 2D and the momenta p and p' thus move 
on either parallel or perpendicular Fermi surface sec- 
tions. In the parallel case, defining p± = (p x ±p y )/2, 
p' ± = (p' x ±p y )/2, we notice that either k + = p + + p' + 
or fc_ = p_ + p'_ vanishes, and the sum in Eq. ((?]) 
(taking, e.g., k + — 0) therefore reduces to x 2£) (a;) ~ 
(4tiV(^))- 1 E q /(e3(q))/(cosg + [cos(^ + fc_)-cos 9 _]). 
Near half filling 8 = A — fij/t is small and the momenta 
effectively extend to < it/2 and \q + \ < it/2 ~ 5, and 
the remaining integrals as well as the averaging of fc_ 
over the square Fermi surface can be calculated. Sum- 
ming over q + gives a logarithmic singularity proportional 
to the level density so that x 2D (At 7 //j) ~ tt/A as fij — >_At. 

Including the induced interactions U e ^ = U\i + Umd 
leads to a gaps simply replacing U\2 by U e ^ in Eqs. (3- 
6). When the induced interactions are small compared 
to the on-site interactions, which is the case for weak 
interactions except near half filling in 2D, we obtain 



A = Aq exp 



3=3 
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where Ao is the uncorrected gap as e.g. given in Eqs. (3- 
5). In a two-component systems the suppression factor 
is consistent with previous estimates of many-body 
effects dUl. 

In an ultracold 6 Li, 40 K or heteronuclear gases the 
many hyperfine states generally have different scattering 
lengths depending on magnetic field. Varying the lat- 



tice depth Vq and wavenumber k allows furthe r tun ing of 
the resulting onsite couplings Uij = £^a,y/c*\/8/7r£ 3 and 
hopping t = E R (2/y/TT)£ 3 e~~ 2 Z 2 , where E R = k 2 /2m is 
the recoil energy and £ = (Vo/Er) 1 / 4 , and subsequently 
also the effective interactions included in Eq. fTTj) , 

In gases with ultracold 173 Yb the scattering len gths 
between the v = 6 nuclear spin states are the same. []J| 
If we assume equal population of all states we find from 
Eq. ifTTj) that the induced interactions change the gap at 
low densities by a factor 



A 2D = A 



2D 



(12) 



in 2D. In 3D A 3D = A|] D (4e) ly / 3 - 1 , since X 3D (^i,^i) = 
ln(4e)/3. The induced interactions from multi- 
components j > 3 enhance superfluidity in such gases 
as well as in symmetric nuclear matter (y = A) and color 
superconductivity in quark matter (8 color, 2 spin and 
2-3 flavor states) [Hj]. 

The above results apply to weakly attractive systems. 
For strong attractions near the unitary limit (BCS-BEC 
crossover) the induced interaction may differ. In a recent 
calculation[25] based on the Nozieres & Schmitt-Rink ap- 
proach the induced interactions reduce pairing by only 
twenty percent whereas experiments Q indicate a reduc- 
tion by a factor ^2 — 3 compatible with an extrapolation 
from the dilute limit including induced interactions [2l| . 

Repulsive on-site interactions U\2 > inhibit s- 
wave pairing unless a longer range attraction is added 
such as a nearest neighbor interaction VJ2(ij) n i n j- 
Super-exchange in the Hubbard model generates a sim- 
ilar nearest neighbor (spin-spin) interaction with cou- 
pling J = At 2 /U\2, which is believed to be responsible 
for high temperature superconductivity [T3|. In the limit 
U\2 — > oo extended s-wave pairing can occur but req uires 
a strong nearest neighbor attraction V < —TT 2 t/2 [lill^. 
whereas d-wave pairing occurs naturally in 2D for even 
weak nearest neighbor attraction. The d-wave mean field 
gap equation is (see, e.g., HF 



1 



V 



A ^ 2E n 
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where now Eq = 



2 + A 2 r] 2 and ?y q = 2[cosq 2; -cosg y ]. 

At low filling the d-wave gap can be calculated within 
mean field to leading orders in the density [171 ] 



A d = 



t 



: exp 



+ c + cin + c 2 n' 



(14) 



The higher order corrections in density are: Co = A/it — 
1 ~ 0.27, ci = tt/2 - 1 ~ 0.57 and c 2 ~ 0.09. 

At half filling we can calculate the d-wave gap within 
mean field as above if correlations can be ignored, which 
requires that the on-site interaction is small. To leading 
logarithmic orders the d-wave pairing gap is 



;t exp 



1 . 



(15) 
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It is a coincidence that the prefactor 8/e 2 is the same 
as in Eq. ([3]). Correlations are, however, expected to 
suppress the d-wave mean field gap of Eq. lfl5|) near half 
filling when the on-site repulsion is strong. 

The influence of induced interactions on anisotropic 
pairing may generally be calculated from the an- 
gular dependence of the screening function xi = 
(x(p/2fei?) cos(Z0)) for 2D multipoles I = 0, 1, 2, .., where 
9 is the angle between pi and P2 (see [13] for the 3D case). 
For a two-component 3D system at low density we find for 
the p-wave Xi D (A i i J Mi) = log(e/4) 1 / 5 ~ —0.08, implying 
that the p-wave induced interactions are attractive and 
enhance p-wave pairing by a factor (4/e) 1 / 5 ~ 1.08. Al- 
though small, it may be important in low density neutron 
matter [13] . In 2D at low density, however, the screen- 
ing function was constant for two components as well as 
multi-components whenever fij > fi\ and since the Fermi 
surface is circular symmetric all higher multipoles than 
the s-wave vanish. Consequently, d-wave pairing is unaf- 
fected by induced interactions in this limit and therefore 
the attractive super-exchange interactions in the spin- 
exchange channel are not suppressed by induced interac- 
tions. However, adding multi-components or approach- 
ing half filling the screening function is no longer sym- 
metric and one can enhance or suppress d-wave super- 
fluidity depending on the sign of the interactions. Near 
half filling the induced interactions are generally large as 



for s-wave and therefore 2D d-wave superfluidity range 
from strongly suppressed for a two-component system to 
possibly strongly enhanced in a multi-component system. 

The critical temperature T c can in 3D be calculated 
from the gap equation ((2} whereas in 2D a Berezinskii- 
Kosterlitz-Thouless transition occurs [29( at a lower tem- 
perature Tbkt- At low density the transition tempera- 
ture is in 3D generally proportional to the gap with the 
same prefactor T c /A = e E /n ~ 0.567.. (E = 0.577.. 
is Euler's constant) for both s- and d-wave superfluidity 
with and without induced interactions [2Sfl. In 2D a mean 
field critical temperature Tmf > Tbkt [23] , can be cal- 
culated from the gap equation as above and we find that 
the same ratio applies even near half filling where loga- 
rithmic singularities appear. This implies that induced 
interactions change the mean field critical temperatures 
by the same factor as the gap of Eq. (fTT]) . 

In summary, ultracold Fermi atoms in traps and op- 
tical lattices provide important systems where we can 
study and test superfluidity and many-body effects in 
detail for multi-component systems. 2D and 3D opti- 
cal lattices have different pairing properties as calculated 
above for gaps, critical temperatures and induced inter- 
actions. The latter can have strong effects on s- and 
d-wave pairing especially near half filling, which affects 
the phases and shell structures near the Mott insulator 
plateaus in traps and confined optical lattices. 
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